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The transdermal drug delivery is the delivery of drugs by using dermal layers. In medical 
sciences, it is an alternative to injections and oral drug delivery system. In this paper a 
mathematical model has been constructed to study the drugs concentration in the different layers 
of the skin through transdermal drug delivery system. The finite element method with linear 
shape functions has been used to obtain the solution of governing one-dimensional partial 
differential equation for unsteady state case. It is assumed that the drug is administered 


internally through reservoir. The dermal region under consideration is divided into five layers 


for the analysis. 
Introduction 


The suitable administration route for strong and low 
molecular weight drugs is transdermal drug delivery. It is 
an substitute to tablets and injections. This delivery 
system is mainly concerned with the delivery device and 
anatomy of dermal region, which consists of uppermost 
stratum corneum and underlying layers of stratum 
germinativum, dermis and subcutaneous tissue [1]. The 
transdermal drug delivery (TDD) is in suitable format is 
supplied externally either through a reservoir in contact 
with the outermost layer or through periodic application. 
The drug and the delivery system are designed in such a 
way that the drug reaches the targeted area in prescribed 
concentration. 

Many researchers have studied drug absorption by two 
methods i.e. experimental and theoretical methods. These 
methods are considered for drug transport generally 
through the transdermal device and the skin. This type of 
transport of drug follows the Fick’s law of diffusion by 
means of a simple homogeneous membrane. Kalia and 
Guy [2] have developed mathematical models based on 
the solution of Fick’s Second law. These models 
described the drug release from various drug delivery 
systems and formulations. These models were based on 
only the different ways of drug transport and analytical 
solutions were found for the drug concentration that 
released from the transdermal device. They considered 
the dermal regions as boundary conditions only. 
“References [3, 4] were the first those have established 
the mathematical drug transport model for one- 
dimensional and single component drug penetration. The 
solutions were obtained by using finite difference 
methods. They have also considered the two pathways i.e. 


intercellular and transcellular for drug penetration in the 
skin. In these models the skin was divided into two sub 
layers as stratum corneum and epidermal layers. The drug 
diffusion model with partition of compounds in biological 
tissue was solved by Missel [5] by using finite element 
method. Chandrasekaren, Micharls, Compbell and J.E. 
Show [6] have examined the tendency of skin during the 
drug permeation process. They obtained the solution of 
transdermal diffusion of scopolamine for single 
component. A notable research work in the direction of 
heat distribution in human dermal region had carried out 
by Saxena and their co-researchers [7-9]. The finite 
element method has been used in most of the cases to 
study the distribution of temperature in the different 
layers of the skin. Gurung [10] studied the time 
dependent temperature distribution in human skin. The 
cold effect in human dermal layers is studied by Khandy 
and Saxena [11]. Using the same approach Khandy [12] 
computed the drug absorption at the nodal points of three 
dermal layers by taking the absorption as a function of 
field concentration. Sharma and Saxena [13-15] have 
established the drug distribution models in transdermal 
drug delivery systems by making use of FEM with 
different shape functions in epidermis, dermis and 
subcutaneous tissues. Also, they have studied the 
transdermal drug flow in radial and angular directions of 
human limbs. 

In fact, the transport process of drugs in transdermal drug 
delivery systems involves the multi layers of the skin. In 
these systems the drug present in reservoir at outermost 
layer stratum cornenum, reaches to targeted region via 
different layers of the skin. Following the above 
approaches, this paper develops a five layered drug 
delivery model with the application of one dimensional 
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variational finite element method to account for effect of 
the interfacial kinetics on the drug delivery process. This 
model is based on the diffusion for transport across a 
homogeneous medium. The present model involves the 
change in material properties, initial conditions and 
boundary conditions in only one dimension. There is no 
*first pass metabolism' effect occurs in TDD therefore we 
ignore the presence of metabolic reactions in the skin. 


Experimental 


Material and Method 
The transport of drugs in the skin is governed by the 
following differential equation taken by Crank [16] 


THOME: 
t 
X X (1) 


where D ,C and ‘ denote the diffusivity, concentration of 
drug and time respectively. Also, 7 represents absorption 
rate of drug by tissues, P is the drug intake rate by blood, 
both are different in each layer, and * is the space 
variable denoting depth below the skin surface. 

The model under consideration, involves layers of the 
skin Figure 1 viz, stratum corneum, stratum 
germinativum, papillary, reticular and subcutaneous 


region with outer and internal boundaries equal to “0 and 


as Tt is also supposed that the applied drug has 


partitioned the uppermost layer stratum corneum and 
entered into the epidermis. The drug concentration is 
defined as a function of one space variable ~ and varies 
linearly in each layer. 
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Figure. 1 General structure of the skin microscopic cross 
section) [17]. 


The concentration of drug varies linearly in each layer 
and is defined as a function of one space variable * . The 
parameters like diffusivity and absorption coefficients are 
assigned constant but different values in respective layers. 
The physiological barriers to the drug distribution at the 
interfaces are also considered. This partitioning of region 


and assignment of values to various quantities are based 
on the biological properties and on the geometrical details 
of the different skin layers. 


In each layer field drug concentration CO G-1234,5) 


is approximated by a linear shape function in Х- 


direction only as 
(i) 


C =а; +8; х (2) 


7i and Bi co-efficient can be written as are 
С,,(0=1,2,3,4) ag 


where the 


given in terms of nodal concentration 


1 a, b 1 C 1 (3) 


f Й] Й E ka 
l a; В; С, (i = 2,3,4,5) (4) 


The drug concentration at the outer surface is known as 


Co and it is considered as the periodic function 


Co =A(l+sin uzt) (5) 
with initial conditions : 
Colao, 0) =A (6) 
С(а,0/-0 (-1,2,3,4,5) (7) 
where ап € 760795235. ае the nodal 
concentrations. "Equation" (5) indicates the periodic 


administration of the drug at the skin surface. Also, ^ 


and / are known as control parameters which determine 
the periodicity and quantity of drug applied. "Equation 
(6)" and “(7)” show that initially the drug is in the 
reservoir, which is situated at the outer surface. At the 
innermost boundary the concentration 
С; =0 

(8) 

It is assumed that the applied drug is targeted to dermis 
only and it has negligible concentration beyond the 
subcutaneous region. This effect can be shown by “(8)”. 


. . . . (і) 1 ЕЕ 
The concentration distributions C",6-212,3,4,5) are 


more or less matching at the interfaces due to mild 
barriers. However, the flux is generally incessant in many 
cases. For this reason, we have 


(0) _ (Hl) - jem 
CY’ =p,C at «^ 4» (i = 1,2,3,4) (9) 


where 2776 are the skin partition coefficients for the drug 


at the respective interfaces. In most cases values of 7 


(= 1,2,3,4) are close to 1. 


Solution 

Many researchers solved such type of problem. The 
boundary value problem for various situations is solved 
by numerical methods because of non-availability of 
exact and analytical solutions. The finite element method 
is one of the most dependable methods in numerical 
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analysis. The usefulness of the method just because of its 
position dependent properties of parameters and their 
flexibility. As the problem involves irregular geometries, 
the method is suitable to understand the feasible 
diffusion. The variational integral 


І-| F(C,C'x)dx 
J (10) 
in optimum form is equivalent to Euler-Lagrange 
differential equation given by Myer [18] 
OF alee iz OC. 
OC dx\cCc' Ox (1 D 
On comparing “(1)” with Euler-Lagrange “(11)", 
obtain the variational form 


i-[ precy Brc 1967 dx 
o| 2 \ Ox 2 Ot 


where the region 'X? is divided into five sub regions. 
1; AG = 1, 2, 3, 4, 5) 


(12) 


The integrals are evaluated for 


eachsubregion and are given as 
DB (T.+B, 

I; M" B, + ( Li l ) 

2 а, 


і 


IL 4a) 4 Pg, + B.) (a; a7) 


ET (а; – a; 4) + За, В,(а? -ар) B; (a; =a | 
(13) 


On using the values of ^; and Bir we get li in terms of 


nodal concentrations Ci, as follows: 


D; (Ci -C? 


l= (аа) + (T, + By (Coa — аус) 


bod 
--- 3(C,a, - aC, 
6(a, —a,) dt 


+3(C, - C,)(G, AC, Xa, -а,) | 


1 
25 + By (ay + ag) (Cy — Co) 


253 Be uud 
E (C, - Со) (af +a + аад) 


(a, — ag) (14) 
е ү 4a; 
F pew? + (T, + B;)( Cat Са, |) 
Ха,-а;|) Pia 


sho, і Say. +a) 
2 p 


2 
P d gest os 
6(a; —aj 4) dt Pia 
2 
si Xa, +a, ) | С, ла, Сам fe = es) 
(a; = d, |) Pia Pia 


(а, +a? +4,4,, ) (с с, | 
(а,-а,|) "Әз А (i = 2,3,4,5) (15) 


Thus integrals I; ( =1, 2,3, 4,5 


+ 


) are assembled to obtain 


i=l (16) 
Now ‘Z’ in “(16)” is extermized with respect to each Ci 


Č =1, 2,3, 4) to given the following system of differentia 


equations: 


4 
> DAD e c 
j-0 (65523,4) (17) 
where PU and Ф/  (/[05234j-01234, are 
constants depending upon various physical апа 
physiological parameters which are listed in Appendix. 
On taking Laplace transform of “(17)” we get 


In, із із 0 (|с, т; 


0 1; із lh C; m3 


(18) 


C; G-L2.5.4) ае the values of Laplace 


where 


. ph 
transform of ©. The constant coefficients — 


(i=1, 2,3,4, j=1,2,3,4) ang "i 6752,54) are listed 
in Appendix. 

On solving “(18)” we get the nodal concentrations as 
= (F 5 +F 3s? + Pas Е; 4) 


ему e Nus) E Nas? Nas ENS) 


Rs OR gs Ra +R;45+R;5 


s? 3a) (М184 +Ns° +N38° +N4s +N; 
(і -1, 2, 3, 4) (19) 


L(C;)- C; = [ceca Ве(ѕ) > 0, ((212,3,4) 
where 0 


А F.. R.,, 
Неге “9” is the Laplace parameter and all '/ and '' ( 


і=1,2,3,4 and. J=1234 are the constants depend 


upon the previously defined parameters. Applying the 
inverse Laplace transform on “(19)” and by making the 
use of Haviside Expansion theorem [19], we get all the 
nodal concentrations as follows: 


3| Uu тақыр TU SET. ) le 


F. 
C =+ 
(Nr +4N,r +3N,r +2М№к + №.) 


Ns 


S| 


(i 21,2, 3, 4) 


(Rar? + т + Rar; + Ra) 
eh 
(Gyr? € Gor? Gr) +Gyr? бық +6) 


where 7176 are the roots of equation 
N,S* + NS? + NS? - NS +N; =0 
N, (і-1,2,3,4,5) 


where are the constant coefficients. 
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Results and Discussion 


For obtaining the drug concentration profile in the 
considered layers, we employ the following physical and 
physiological parameters which are practicable and fall 
within the acceptable ranges. We can allocate any value 
to the thickness of in-vivo tissues depending on the 
sample of the skin under study. Accordingly the set of 
values considered here are 

а =2.0 ит, а,-4.0шт, 


ag = 0.0 um., аз =8.0 um 


a4 =12.0umm а =16.0 um 


p;-Lh (і-1,2,3,4) 


4-100. и= 0.01. 


The mathematical calculations have been carried out for 
different values of applied drug concentration, drug rate 
of absorptions and diffusivity. The terms T and В are 
taken equal to zero in upper most layer because this layer 
is free from active cells and blood vessels and have 


(i) 
different values in other layers. The graphs for с, 


( — 1, 2, 3,4, 5) versus thickness of the skin layer and time 


have been drawn using MATLAB software. 

The present work is an attempt to estimate the drug 
concentration in human dermal layers. Firstly we have 
approximated the drug concentration at the nodal points 
of the layers of the skin with respect to time. It has been 
observed from the curves given in figure 2, 3 and 4, that 
the steepness of the curves decreases as the drug reaches 
to internal boundary. This effect is due to the boundary 
conditions where concentration decreases with raise in 
distance from the surface. The effect of change in the 
values of absorption coefficients and diffusivity can also 
be observed from these curves. The increase in absorption 
decreases the concentration. This effect is shown in figure 
2 and 3. 
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Also, the increase in diffusion increases the concentration 
is shown in figure 3 and 4. Figure 5 to 7 are drawn 
between field concentration and depth for different values 
of time. The slop of these curves decreases at the 
interfaces. This result is due to the variation in properties 
of each sub section. The curves are linear within each 
element due to the linear shape functions, which are 
supposed in each element. Three-dimensional “Figure 8 
to 10" are also drawn to show the variations in drug 
concentration with respect to time and thickness of 
dermal layers at the same time. 
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Conclusion 


The solutions provided in this paper give better 
understanding of mechanisms of drug delivery in the five 
layers of the human skin. The results in this paper are 
stronger than the results obtained by Khandy and Rafiq. 
They have obtained the concentration only in three layers 
while as in our model, the drug concentration profiles has 
been studied for five layers as the drug flow in the skin is 
multi layered observable fact. Also, they have studied the 
drug distribution in-vivo tissues by applying the drug 
amount 5 grams while as the present work shows the drug 
distribution for different amount of applied drug 1.e. 100 
grams, 160 grams and 200grams. Hence our work gives 
modified results followed by the Khandy and Rafiq. 
Moreover, our work shows the linear variation in each 
element. It is a justified approximation in view of small 
thickness of the layers. This assumption can always be 
enhanced by increasing the number of layers or by taking 
higher degree approximating polynomials. The developed 
model helps us in calculating the amount of drug 
concentration in different layers of the skin with the effect 
of change in diffusion constant and absorption 
coefficients. The model may be very helpful in 
formulation the effective transdermal drug delivery 
systems in various situation. 


Appendix 
_ D; = 
Pio = 0, ат pasa) | En pi370-7 ра, Фо 70, 
di? = (a5 –а)/6р 913 20 7 44, P20 =0, 
Р, 7 
2 Рэ(аҙ-а;) ius P24 = 0, 40 =0= 44, 
(а; —a5) 
HEN бр; ES P30 707 P31, 
P34 = 1 = P43 
p3(a4 — 43) : D34 — —D,/ p (a4 —43) = P43, 
934 = 2a = 443 
dij =0= 431, брз Рад 707 ра = ра, 


440 = 0= 441 7 942, l3 = 1 754,251) ті 71570, 
lj = Ру +595 (і-1,2,3,4,) =1,2,3,4), 
_ (T; t Bj) (a; - aj 1) b (Т, Bia) (aj; - aj) 


h; 
2 2 pi 
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ne 13 " hy 
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Е DA " 
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12 = ЕЕ 
2 2p, 
DAx Di Dj 
£7— ux рі = toz - 
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qi = (a; —a; 4) " (Gist 80 - 
3 3 pi (i =1, 2, 3, 4) 
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